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Resumen 

La colision de una onda de choque piano paralela con una nube piano paralela de densidad 
uniforme es analizada para el caso en el que campos magneticos y perdidas por radiacion no son 
consideradas. Se discuten soluciones anahticas generales para el caso en el que la densidad de la 
nube es mucho mayor que la del gas que le rodea. Este problema generaliza uno de los problemas 
clasicos en el estudio de la dinamica de gases: la colision entre una onda de choque y una pared 
solida. 

Abstract 

The collision of a plane parallel shock wave with a plane parallel cloud of uniform density is 
analysed for the case in which magnetic fields and radiative losses are not considered. General 
analytic solutions are discussed for the case in which the density of the cloud is greater than that 
of the surrounding environment. This problem generalises one of the classical problems in gas 
dynamics: the collision between a shock wave and a solid wall. 

PACS:47.40.N 

§1 Introduction 

The problem of the collision of a shock wave with a cloud has been intensely investigated in the 



past by several authors (see for example Klein et ah 1994 and references therein). The simplest 



assumption to make is to consider a cloud for which gravitational effects are not considered, magnetic 
fields are non-important and radiative losses are negligible. The fact that gravity is not taken 
into account, makes it possible to consider the density of the cloud as uniform. The complete 3D 
hydrodynamical problem is extremely complicated, even under the simplifications mentioned above. 
However, numerical simulations have been done for this case which ultimately give rise to instabilities 



causing a complete disruption of the cloud ([Klein et al., 1994|) . 



This article describes how the solution of the one dimensional problem can be obtained. It has 
been argued in the past that at least for the case in which the density contrast is high, i.e. the ratio 
of the cloud's density to that of the external environment is high, the problem has to be very similar 



to the one found in the problem of a collision of a plane parallel shock with with a solid wall (Spitzer 
198^ , |McKee 1988; ). 



Many Astrophysical phenomena give rise to collisions between a shock wave and a cloud. For 
example, when a supernova explosion occurs, the intense ejection of energy from the supernova into 
the interstellar medium produces a spherical shock wave which expands into the interstellar medium. 
Several examples exist for which collisions of this expanding shock have been observed to interact 



with clouds embedded in the interstellar medium. This interaction is very important, since it seems to 
induce, under not very well known circumstances, gravitational collapse and star formation ( [Herbst 



an i Assousa, 1979| ). Another scenario is presented by the expansion of jets around active galactic 
nuclei. A pair of jets expand in opposite directions from the nuclei of the galaxy creating a bow 
shock which interacts with the intergalactic medium. It is the interaction of this expanding bow 
shock with clouds or galaxies embedded in cluster of galaxies that provides a mechanism in which 
shock-cloud interactions take place. It seems that this interaction is able to induce star formation 
very efficiently^. 

Having all this considerations in mind, the present paper aims to give a simple way of solving 
a particular case of the whole problem. This article provides an analytic description of the one 
dimensional problem of a collision between a plane parallel shock with a plane parallel "cloud" 
bounded by two tangential discontinuities. It is assumed that the specific volume in the cloud is a 
quantity of the first order, in other words solutions are given for the case in which the density of the 
cloud is much greater than that of the surrounding environment. 



§2 General Description of the problem 

Consider two plane parallel infinite tangential discontinuities. The cloud, or internal region to the 
tangential discontinuities has uniform pressure p c and density p c . The environment, or external region 
to the cloud has also uniform values of pressure p\ and density p\ respectively. A plane parallel shock 
wave is travelling in the positive x direction and eventually will collide with the left boundary of the 
cloud at time t = to < 0. For simplicity we assume from now on that the density of the cloud is greater 
than that of the environment. By knowing the pressure P2 and density pi behind the shock wave, it 
is possible to solve the hydrodynamical problem thus defined. 



The problem of the collision of a shock wave and a tangential discontinuity is well known (Landau 



and Lifshitz, 1987 ). Since at the instantaneous time of collision the values of, say, the density in 
front and behind the shock are p c and p2 respectively, the standard jump conditions for a shock no 
longer hold. A discontinuity in the initial conditions (first initial discontinuity) occurs. 

When a discontinuity in the initial conditions occurs, the values of the hydrodynamical quantities 
need not to have any relation at all between them at the surface of discontinuity. However, certain 
relations need to be valid in the gas if stable surfaces of discontinuity are to be created. For instance, 
the Rankine-Hugoniot relations have to be valid in a shock wave. What happens is that this initial 
discontinuity splits into several discontinuities, which can be of one of the three possible types: shock 
wave, tangential discontinuity or weak discontinuity. This newly formed discontinuities move apart 
from each other with respect to the plane of formation of the initial discontinuity. 

Very general arguments show that only one shock wave or a pair of weak discontinuities bounding 
a rarefaction wave can move in opposite directions with respect to the point in which the initial 
discontinuity took place. For, if two shock waves move in the same direction, the shock at the front 
would have to move, relative to the gas behind it, with a velocity less than that of sound. However, 
the shock behind must move with a velocity greater than that of sound with respect to the same gas. 
In other words, the leading shock will be overtaken by the one behind. For exactly the same reason 
a shock and a rarefaction wave can not move in the same direction, and this is due to the fact that 
weak discontinuities move at the velocity of sound relative to the gas they move through. Finally, 
two rarefaction waves moving in the same direction can not become separated, since the velocities of 
their boundaries with respect to the gas they move through is the same. 

Boundary conditions demand that a tangential discontinuity must remain in the point where the 
initial discontinuity took place. This follows from the fact that the discontinuities formed as a result 
of the initial discontinuity must be such that they are able to take the gas from a given state at one 
side of the initial discontinuity to another state in the opposite side. The state of the gas in any one 



See for example the Hubble Space Telescope WWW site at tittp: / / oposite.stsci.edu /pubinfo /pr /1995/ 30.htm] 



dimensional problem in hydrodynamics is given by three parameters (say the pressure, the density 
and the velocity of the gas). A shock wave however, is represented by only one parameter as it seen 
from the shock adiabatic relation (Hugoniot adiabatic) for a polytropic gas: 



= (7 + 1)P/ + (7 ~ l)Pb (2 n 

V f (7-l)P/ + (7 + l)Pb' 

where p and V stand for pressure and specific volumes respectively, 7 is the polytropic index of 
the gas and the subscripts / and b label the flow in front of and behind the shock. For a given 
thermodynamic state of the gas ( i.e. for given pf and Vf) the shock wave is determined completely 
since, for instance, pb would depend only on V& according to the shock adiabatic relation. On the 
other hand, a rarefaction wave is also described by a single parameter. This is seen from the equations 
which describe the gas inside a rarefaction wave which moves to the left with respect to gas at rest 
beyond its right boundary ( Landau and Lifshitz, 1987|) : 



cr = c 4 + -(j c -1)w r , (2.2) 

w R = 2 —rU + j). (2.5) 

7 C + 1 V tJ 

where C4 and cr represent the sound speed behind and inside the rarefaction wave respectively. The 
magnitude of the velocity of the flow inside the rarefaction wave is wr in that system of reference. 
The quantities p^ and pr are the pressures behind and inside the rarefaction wave respectively. The 
corresponding values of the density in the regions just mentioned are p± and pr. 

With only two parameters at hand, it is not possible to give a description of the thermodynamic 
state of the gas. It is the tangential discontinuity, which remains in the place where the initial 
discontinuity was produced, that accounts for the third parameter needed to describe the state of 
the fluid. 

When a shock wave hits a tangential discontinuity, a rarefaction wave can not be transmitted to 
the other side of the gas bounded by the tangential discontinuity. For, if there would be a transmitted 
rarefaction wave to the other side of the tangential discontinuity, the only possible way the boundary 
conditions could be satisfied is if a rarefaction wave is reflected back to the gas. In other words, 
two rarefaction waves separate from each other in opposite directions with respect to the tangential 
discontinuity that is left after the interaction. In order to show that this is not possible, consider a 
shock wave travelling in the positive x direction, which compresses gas 1 into gas 2 and collides with 
a tangential discontinuity. After the interaction two rarefaction waves separate from each other and 
a tangential discontinuity remains between them. In the system of reference where the tangential 
discontinuity is at rest, the velocity of gas 2 is v 2 = — P 2 \/—dpdV, where p% is the pressure of gas 
3 surrounding the tangential discontinuity. Accordingly, the velocity of gas 1 in the same system of 
reference is v\ = — \J — dpdV . Since the product — dpdV is a monotonically increasing function of 
the pressure and < pj, < p% then: 

rP2 rpi rp2 

- / V-dPdV <v 1 -v 2 < / V-dPdV - / V-dPdV. 

JO JO Jpi 

The difference in velocities v% — v 2 has the same value in any systems of reference and so, it follows 
that v\ <t>2, in particular on a system of reference with the incident shock at rest. However, for the 



incident shock to exist, it is necessary that v\ >V2, so two rarefaction waves can not be formed as a 
result of the interaction. 

So far, it has been shown that after the collision between the shock and the boundary of the 
cloud, a first initial discontinuity is formed. This situation can not occur in nature in any manner 
and the shock splits into a shock which penetrates the cloud and either one of a shock, or a rarefaction 
wave (bounded by two weak discontinuities) is reflected from the point of collision. With respect 
to the point of formation of the initial discontinuity, boundary conditions demand that a tangential 
discontinuity must reside in the region separating the discontinuities previously formed. 

In a shock wave, the velocities (v) in front and behind the shock are related to one another by 
their difference: 



Vf-v b = y/(p b ~ PfW f - H), (2.6) 

where the subscripts / and b label the flow of the gas in front and behind the shock wave. 

If after the first initial discontinuity two shock waves separate with respect to the point of collision, 
then according to eq. (p.6[) the velocities of their front flows are given by v c = — \J (p^ — pi)(V c — Vy) 
and V2 = v/ (p3 — P2) ( V2 — V3) , where the regions 3 and 3' bound the tangential discontinuity which 
is at rest in this particular system of reference (see top and middle panels of fig. ([!])). Due to the fact 
that pz > P2 and because the difference V2 — v c is a monotonically increasing function of the pressure 
P3, then: 

v 2 -v c > (p 2 -piW{2V c / [( 7c - I) Pi + (7c + 1)P2]} , 
according to the shock adiabatic relation. Since V2 — v c is given by eq.( |2.6D , then: 



YL > Ye (27) 

(7 - 1) + (7 + IWpi (7c - 1) + (7c + IWpi ' 

where 7 and 7 C represent the polytropic indexes of the environment and the cloud respectively. V\ and 
V c are the specific volumes on the corresponding regions. In other words, a necessary and sufficient 
condition for having a reflected shock from the boundary of the two media, under the assumption of 



initial pressure equilibrium between the cloud and the environment, is given by eq.(2.7). Since for 
the problem in question V\>V C and the polytropic indexes are of the same order of magnitude, a 
reflected shock is produced. 

In the same form, at time t = when the transmitted shock reaches the right tangential discon- 
tinuity located at x = 0, another (second) initial discontinuity must occur. In this case, we must 
invert the inequality in eq,(|2.7j), change 7 by 7 C and P2 by p^, where p% is the pressure behind the 
shocks produced by the first initial discontinuity. Again, using the same argument for the poly- 
tropic indexes, it follows that after this interaction a weak discontinuity bounded by two rarefaction 
waves must be reflected from the boundary between the two media. As a result of the interaction, 
once again, boundary conditions demand that a tangential discontinuity remains between the newly 
formed discontinuities. 

This situation continues until the rarefaction wave and the left tangential discontinuity of the 
cloud collide at time t = r > 0. At this point, two rarefaction waves separating from each other 
from the point of formation will be produced once a stationary situation is reached, and a tangential 
discontinuity will be separating the newly formed discontinuities. One can continue in a somewhat 
indefinite manner with the solution but, for the sake of simplicity the calculations are stopped at this 
point. Fig.(|]]) shows a schematic description of the solution described above in a system of reference 
such that the tangential discontinuities which are left as a result of the different interactions are at 
rest. The numbers in the figure label different regions in the flow. 
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Figure 1: An incoming shock travelling to the right (top panel) hits a tangential discontinuity at time to <0. This 
produces two shocks moving in opposite directions with respect to the place of formation (middle panel). When the 
transmitted shock into the cloud (region C) collides with its right boundary a reflected rarefaction wave (region R) 
bounded by two tangential discontinuities and a shock transmitted to the external medium (lower panel) are formed. 
Arrows represent direction of different boundaries, or the flow itself. The numbers in the figure label different regions 
of the flow. Dashed lines represent shocks, dash-dot are weak discontinuities and continuous ones are tangential 
discontinuities. The system of reference is chosen such that the tangential discontinuities which are left as a result of 
the collisions are always at rest. 

§3 First initial discontinuity 

According to fig.(|l]), after the first initial discontinuity the absolute values of the velocities (v) of the 
flow are related by: 

V 2 + V c = V 2i . (3.1) 

With the aid of eq. (|2.6D , the velocities of eq (|3.1[) are given by: 



vl = (p2-pi)(V 1 -V 2 ), (3.2) 
v 2 c = (P3-Pi)(V c -V 3 ,), (3.3) 

(3.4) 
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Inserting eqs. (|3.2| )- (l3.4j) into eq. (|3.l|) and substituting for the specific volumes from eq. (|2.l|) , one 
ends with a relation which relates the pressure p 3 as a function of p 2 , pi and the polytropic indexes 
in an algebraic linear form. Straightforward manipulations show that the resulting equation does 
not have an easy analytic solution, even for the particular cases in which a strong or weak incident 
shock collides with the cloud. 

In order to find a set of analytic solutions, let us first describe a particular solution to the problem. 
If we consider a cloud with an initial infinite density -a solid wall, then eq.(3.1) takes the form v 2 = v 2 i, 
and a "zeroth order" solution is found ( [Landau and Lifshitz, 1987| ): 

P3o = (37 - l)p 2 ~ (7 ~ , 3 ^ 

P2 (7 - 1)P2 + (7 + ' 

where p 3o is the value of the pressure behind the reflected and transmitted shocks for the case in 



which the cloud has specific volume V c = 0. For this particular case, eq^^) determines ps as a 
function of p\ and p 2 , which are initial conditions to the problem. Due to the fact that the gas is 
polytropic, this relation is the required solution to the problem. 
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In order to get a solution more adequate to the general case, we can approximate the whole 
solution under the assumption that V c is a quantity of the first order, so: 



P3 =P3 +P3> 
V 3 = V 3o +Vf, 



(3.6) 
(3.7) 
(3.8) 



where the quantities with a star are of the first order and the subscript represents the values at 
zeroth order approximation. Substitution of eqs.( |3.6j )-( |3.8| ) into eqs.(^)-( |3.4| ) gives: 



vl 



Ac 



-V 3 *(p 3o -p 2 )+p* 3 (V 2 -V 3o ), 

c = (P3 -Pl)(V c 



From the shock adiabatic relation, eq.(2.1), and eqs. (|3.6i )- (|3.8| ) it follows that 

(7 + l)p 2 + (7 - l>3o 
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V 2 
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(7 - 1)P2 + (7 + 1)P3 ' 
(7c + + (7c ~ l)P3o 

Vc (7c - + (7c + 1)P3 
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Substitution of eqs.((3 



v 3 = 

V 2 ' [( 7 _i) P2 + ( 7 + i)p 3o ] 2 - 
-( |3.1C| ) and eq.( |3.13| ) in eq,( |3.1| ) gives the required solution: 



(3.9) 
(3.10) 



(3.11) 
(3.12) 
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where: 
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(3.14) 



The specific volumes V 3o and Kj* are given by eq.(3.11) and eq.( |3.12 ) respectively. For completeness, 
approximations to eq,( |3.14D for the case of a very strong incident shock and that of a weak incident 
shock are given: 



where: 
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(3.15) 
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Figure 2: Variation of the pressure P3 behind the transmitted shock into the cloud as a function of the strength of 
the initial incident shock. The continuous line shows the case for which the cloud is a solid wall with infinite density. 
The dashed curve is the solution at first order approximation for which the cloud's specific volume is a quantity of 
the first order. The acoustic approximation for which the incident shock is weak, at the same order of accuracy, is 
represented by a dot-dashed curve. The perturbed solutions were plotted assuming p c /pi = 100 for polytropic indexes 
7 = 7 C = 5/3, corresponding to a monoatomic gas. 



and C= (p2—pi)/pi "C 1 in the weak limit. Fig.(^) shows a plot of the pressure ps as a function of the 
strength of the incident shock. It is interesting to note that even for very strong incident shocks the 
ratio pz/p2 differs from zero, which follows directly from eq.(|3.5|) and eq.( |3.15| ). This simple means 
that the reflected shock is not strong, no matter the initial conditions chosen. 

There are certain important general relations for which the above results are a consequence 
of. Firstly, by definition the pressure p2 behind the shock is greater than the pressure p\ of the 
environment. Now consider a strong incident shock, then since ps > P2 ^ Pi, it follows that the 
transmitted shock into the cloud is very strong. Also, the reflected shock does not have to compress 
too much the gas behind it to acquire the required equilibrium, so it is not a strong shock. This last 
statement is in agreement with eq.( 3.15| ). In general, for any strength of the incident shock, since 
the inequality P3>P2>Pi holds, continuity demands that the reflected shock can not be strong and, 
more importantly, that the penetrating shock is always stronger than the reflected one. 

Secondly, very general inequalities are satisfied by the velocities V2, v c , v s i as defined in fig. ([[]). 
For instance: 



v sl >v 2 . 



(3.17) 



Indeed, since V3V2/CV3 — V2) >V2 holds, and the left hand side of this inequality is just v s i according 
to mass flux conservation across the reflected shock, the result follows. 

On the other hand, from eq.( |3.3| ) and eq,( |3.4| ), since P2 > Pi it follows that a necessary and 
sufficient condition for V2 > v c to be true is that V2 — V3 > V c — Vy. This last condition is satisfied 
for sufficiently small values of V c . To give an estimate of the smallness of the cloud's specific volume 
needed, note that a necessary and sufficient condition for V2 — V3 > V c — Vy to be valid is: 
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(3.18) 



according to the shock adiabatic relation for the transmitted and reflected shocks. Since ^3 >P2>Pi 
and V2 < V\ it follows that: 



Vi 



> 



Vr 



(7 - 1) + (7 + I)p 3 /Pl (7c - 1) + (7c + IWpi ' 



(3.19) 



which is very similar to eq.( |2.7| ). In the same fashion, under the assumption that the polytropic 
indexes are of the same order of magnitude, eq.( |3.19| ) implies V C <V\, which was an initial assumption. 
Although eq.( 3.19j ) is not sufficient, due to the fact that V c is a first order quantity we can use in 
what follows: 



v 2 > v c . 



(3.20) 



The inequalities in eq. (ETTA and eq. (^20h will prove to be useful later when we choose a more 
suitable reference system to describe the problem in question. 



§4 Second initial discontinuity 

Let us now analyse the situation for which 0<i<r. To begin with let us prove that: 



Wi < V 2 + V c = U 2 , 



(4.1) 



where the velocities w\, V2 and v c are defined in fig.(|l|). Suppose that the inequality in eq.(4.1) is not 
valid, then, by expressing the velocities as function of the specific volumes and pressures by means 



of eq.(2.6) and the fact that p 2 >Pi, P3 >Pa and V41 > V3, it follows that p% > p c ; then as the cloud's 
density grows without limit, so does p^. Necessarily, eq.( |4.1| ) has to be valid for sufficiently small 
values of the cloud's specific volume. It is important to point out that since W2 = \u2~ W\\ =112— wi, 
the gas in region 2 as drawn in fig.(Q) travels in the positive x direction. According to fig.(|l|) flows 
in region 1 and 3 are related by 



(4.2) 



Let us now prove a very general property of the solution. Regions 2 and 3 are related to one 
another by the shock adiabatic relation. Since the gas in regions 3' and 4 obey a polytropic equation 
of state Pz/pa = iyi/y? J l V c it follows that: 
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Vy V c (7 + 1)^2 -(7 -1)V3 
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Now, due to the fact that V3/ < V4 < V\ V3 < V 2 < V\ and 7, 7 C > 1 for a reasonable equation of state, 
this relation can be brought to the form 



P4 1 
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Pi 

P2 



(4.3) 



according to the shock adiabatic relation. This result implies that most of the energy from the 
incoming shock has been injected to the cloud, no matter how strong the initial incident shock is. 
Only a very small amount of this energy is transmitted to the external gas that lies in the other side 
of the cloud. Note that this result is of a very general nature since no assumptions about the initial 
density contrast of the environment were made. 

In order to continue with a solution at first order approximation in V c , note that we have to use 
eqs.([3lf)- (l3.8| ) together with: 

P4=Pi+pt, (4.4) 
V 4 = V:, (4.5) 
V 4 , =Vx + V*,, (4.6) 

where the quantities with a star are of the first order. The velocities w\ and w% can be expressed as 
functions of the specific volumes and pressures by means of eq. (|2.6j ), for which after substitution of 
eqs.Q-g]]) it follows: 

vo\ = -p%V£, (4.7) 

^3 = ~^-r (VlcP3 V£ - y/lcPlVA . (4.8) 

7c J- ^ ' 
The specific volumes behind the transmitted shock and the reflected rarefaction wave are obtained 
from the shock adiabatic relation and the polytropic equation of state for the gas inside the rarefaction 
wave: 

V* = - Vl tiL t (4.9) 
V£ = V£, ° . (4.10) 



By substitution of eqs.( f4.7| )-( [4.10 ) and eq.( |3.4| ) in eq.(fT^) the required solution is found: 



P2 y P2 V\ 
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For completeness the limits for the case of strong and weak incident shocks are given: 



P2 



7 (3 7 - 1) pi V c 
(7 C + 1)( 7 -1)P2^1 



El 

P2 



6C 



(4.12) 
(4.13) 



J. \J 




P1/P2 



Figure 3: Variation of the pressure p4 behind the transmitted shock into the external medium as a function of the 
strength of the incident shock. The continuous line represents the case for which the cloud has infinite density and so 
it does not transmit any shock to the external medium. The dashed curve represents the case for which the cloud's 
specific volume is a quantity of the first order. The dash-dotted curve is the limit for which a strong (or weak) incident 
shock collides with the cloud at the same order of approximation. The perturbed curves were produced under the 
assumption that /9 c /pi = 100 for monoatomic gases. 



with: 

fpi 7 -i \M)/V 

\p 2 37 - 1 J 

It follows from eq.( |4.12| ) that p^ <C P2 as the strength of the incident shock increases without 
limit. This result was given by a very general argument in eq.(fD|). Fig.(||) shows the variation of 
the pressure p^ behind the shock transmitted to the environment as a function of the strength of the 
initial incident shock, after the second initial discontinuity. 



>5 General solution 



Having found values for the pressures p% and p\ as a function of the initial conditions p\, P2, V\ 
and V c , the problem is completely solved. Indeed, using the shock adiabatic relation V2 is known. 
With this, the values of Kj*, V£, V£ and V£, are determined by means of eq-dUg), eq.(fyj), eq.pjJ) 
and eq.( [4.10 ) respectively. The complete values for pressure and specific volumes are obtained thus 
with the aid of eqs. (|3.6|) -(3.8) and eqs.(4.4)-(4.6). The velocities of the flow, as defined in fig.(||), 
are calculated either by mass flux conservation on crossing a shock, or by the formula given for the 
velocity discontinuity in eq. ( |2.6D . The hydrodynamical values of the pressure pn and density pn 
inside the rarefaction wave come from eqs.(|2.2|)-( |2.5| ). 

In order to analyse the variations of the hydrodynamical quantities as a function of position and 
time, let us now describe the problem in a system of reference in which the gas far away to the 
right of the cloud is always at rest, as presented in fig.(Q). Let xu and xt T be the coordinates of the 
left and right tangential discontinuities, x s i and x sr the coordinates of the reflected and transmitted 
shocks produced after the first initial discontinuity, Xsr the position of the transmitted shock after 
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Figure 4: Description of the problem of a collision of a shock with a cloud in a system of reference for which the 
gas far away to the right (at x = oo) is always at rest. Originally a shock is travelling to the right and hits a tangential 
discontinuity (top panel). This produces a discontinuity in the initial conditions so a reflected and transmitted shock 
are produced; the gas in the cloud begins to accelerate (middle panel). Eventually the transmitted shock into the 
cloud collides with its right boundary producing a "reflected" rarefaction wave bounded by two weak discontinuities 
(region R) and a transmitted shock into the external medium (lower panel). In this system of reference every single 
discontinuity produced by means of the interaction move to the right, except for the reflected shock produced after 
the first collision. Arrows represent the direction of motion of various boundaries and direction of flow. Numbers label 
different regions of the flow. Dashed lines represent shock waves, dash-dotted ones weak discontinuities and continuous 
ones tangential discontinuities. 



the second initial discontinuity and x a and Xb the left and right weak discontinuities which bound 
the rarefaction wave. The new velocities are defined by Galilean transformations: 

(5.1) 
(5.2) 
(5.3) 
(5.4) 
(5.5) 



"2 


= v 2 + v c , 


Usl 


= V s l ~ V c , 


u sr 


— ^sr ^C; 




= Wl~ WR, 


V sr 


= W\ + w sr 



The direction of motion of the flow is shown in fig.(||) and it follows from eq.( 3.17| ), eq.( 3.20| ) and 



eq. (|5,2| ) that u s i points to the left in this system of reference. Since, in the same frame, v c and w\ 
point to the right, continuity across a weak discontinuity demands vr to do it in the same way. 

The tangential discontinuities and the shocks produced by the initial discontinuities move with 
constant velocity throughout the gas. This implies that the time at which the first initial discontinuity 
takes place is: 



t = - — , (5.6) 



v 



sr 



where A represents the initial width of the cloud. Hence, the positions of all different discontinuities 
for to < t < are: 



X S r — v sr t^ (5.7) 
x s i = -A - u s i(t - t ), (5.8) 
x u = -A + v c (t-t ). (5.9) 



and for 0<i<r, eqs.( |5.8] )-( |5.9| ) are valid together with 



X a = ~t ^ 7C ^" 1 W 3 + C^j + Wit, (5.10) 

x b = (w 1 -C4)t, (5.11) 
Xsr = v sr t, (5.12) 
Xtr = w±t. (5.13) 

The time r at which the left tangential discontinuity collides with the left boundary of the rarefaction 
wave is given by x t i = x a , and thus: 



TC3 , = Vc t Q + A. 



(5.14) 



Fig. (|5|) shows the variation of the pressure and density as a function of time and position in a system 
of reference in which the gas far away to the right of the cloud is at rest. 

The width of the cloud varies with time, and it follows from eq.(5.9) and eq.( |5.13 ) that this 
variation is given by: 



X(t) = Q(t)wit + A - v c (t - t ), 
where @(t) is the Heaviside step function. This linear relation is plotted in fig.® 



(5.15) 



§6 Summary 

The problem of a collision of a plane parallel shock wave with a high density cloud bounded by two 
plane parallel tangential discontinuities has been discussed. Radiation losses, magnetic fields and 
self gravity of the cloud were neglected. General analytic solutions were found for the simple case in 
which the ratio of the environment's density to that of the cloud's density is a quantity of the first 
order. 

When the shock collides with the boundary of the cloud, a discontinuity in the initial conditions 
is produced. This splits the incoming shock into two shock waves: one which penetrates the cloud 
and one which is reflected back to the external medium. When the transmitted shock into the cloud 
reaches the opposite boundary, another discontinuity in the initial conditions is produced, causing 
the transmission of a shock wave to the external medium and the reflection of a rarefaction wave 
from the point of collision. 
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Figure 6: Variation of the width of the cloud in units of its original size A as a function of the dimensionless quantity 
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